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Thermoelectric coolers (TECs) use the Peltier effect for thermal management of 
electronic devices. They offer high reliability and low noise operation but limited in use 
due to low performance. In the present work, through the use of a genetic algorithm 
(GA), two single-objective optimizations associated with two separate objectives are 
carried out, aiming maximization of cooling capacity and maximization of the 
coefficient of performance (COP) of TEC with space restrictions. Interfacial thermal 
resistance and electrical contact resistance are taken into consideration to obtain a more 
realistic model. This paper presents a new approach to finding appropriate solutions by 
optimally arranging the length of n-type and p-type thermoelectric (TE) elements, the 
cross-sectional area of TE elements, and input electric current. To validate the GA 
predictions, three-dimensional steady-state TEC models are prepared, and finite- 
element simulations are carried out using ANSYS®. Close agreement between the GA 
and ANSYS® has been observed. This study provides a new mathematical optimization 
model that is more realistic and is quite close to the physical construction of TEC 
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modules manufactured by industry. 


1. INTRODUCTION 


The solid-state thermoelectric (TE) technology attract great 
attention of the researchers because of its potential use as 
green energy conversion devices. The Peltier effect of 
thermoelectric technology offers direct conversion of 
electrical energy into temperature difference. Conversely, the 
Seebeck effect of TE technology provides the conversion of 
thermal energy of temperature differential into electric power 
[1]. A thermoelectric cooler (TEC) dissipates the heat and 
removes the hotspots of the electronic devices in an 
environment-friendly manner using the Peltier effect. A TEC 
could be installed easily within a restricted space due to its 
practical manufacturing possibility in small _ sizes. 
Thermoelectric coolers must be appropriately designed and 
manufactured to meet the necessary performance requirements. 
Two essential performance parameters of a TEC are the 
cooling capacity and the coefficient of performance. The 
cooling capacity of thermoelectric coolers ranges from 
milliwatts to watts depending on the requirements. The 
maximum cooling effect or higher COP for a thermoelectric 
cooler can be achieved through upgraded TE materials and 
improved device design. 

The efficiency of TE materials increases with a material 
property known as figure of merit (Z). The term Z is defined 
as «7/RK, where a is the Seebeck coefficient, R is the electrical 
resistance, K is the thermal conductance. With absolute 
temperature (7), the dimensionless figure of merit (Z7) is used 
to characterize TE materials. A higher value of ZT corresponds 
to better cooling performance. Hicks et al. described that the 
value of ZT could be enhanced by reducing the dimensions of 
thermoelectric materials [2, 3]. At room temperature, 
Venkatasubramanian et al. [4] reported a ZT~2.4 for p-type 
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BizTe3/Sb2Te3 superlattice devices. Peak ZT values of different 
TE materials are attainable at different temperatures. Over the 
past two decades, significant progress in maximizing ZT has 
been made in developing thermoelectric materials [5-10]. 

With the significant ongoing efforts to improve TE 
materials, the researchers also focus on designing and 
assembling the TECs. The investigations established that the 
geometric structure of thermoelectric elements affects the 
performance of thermoelectric coolers [11-15]. Huang et al. 
[16] combined a three dimensional TEC model with a 
simplified conjugate-gradient technique. They reported that at 
a fixed temperature difference and fixed current, a substantial 
value of the total area of TE elements with small element 
length can maximize cooling capacity. Yang et al. [17] 
reported that micro-thermoelectric coolers operating in a 
transient regime could provide a better cooling effect. Nain et 
al. [18] reported that a suitable value of dimensionless current 
can enhance the performance of TEC. Pareto-optimal 
solutions were obtained for different settings of temperature 
ratio. Shen et al. [19] reported that a two-segmented TE 
element structure can reduce the joule heating effect from 50% 
to 35% on the cold side. The results showed a remarkable 
118.1% improvement in maximum cooling capacity. Nain et 
al. [20] optimized cooling capacity and COP performance of 
TEC using dimensional structural parameters as design 
variables. The geometrical parameters were optimized to find 
Pareto-optimal solutions. Jeong [21] reported that the COP of 
TEC can be increased by optimal values of current and length 
of thermoelements. Lee [22] proposed a dimensional analysis 
approach to find out the optimal design of TE devices with 
feasible mechanical constraints. Mijangos et al. [23] reported 
a novel design of asymmetrical legs to enhance the 
performance of TE devices. 


Literature reports several studies on performance 
optimization of TEC [13, 15, 24-29]. However, in the current 
study, the two performance parameters, namely, cooling 
capacity and coefficient of performance, are optimized as two 
single-objective optimization problems. So far, the standard 
approach has been to choose either a set of geometric design 
variables or operating design variables. In this paper, a 
combined set of three design variables, electric current, length 
of n-type and p-type TE elements and cross-sectional area of 
TE elements is chosen in both optimization problems. The 
optimization algorithm mathematical model is customized to 
handle the presence of ceramic substrate, copper contacts, 
electric contact resistances at the interface, and heat sink, 
which are essential parts in the fabrication of a TEC module in 
industrial applications. It is a new aspect of modelling TEC. 
The geometry of the thermoelectric element plays a vital role 
in the performance of the thermoelectric cooler. However, 
tight geometric space constraints are found in many 
telecommunications and other scientific applications. The 
TEC is used for cooling electronic devices where space 
restrictions are quite prevalent. Hence, consideration of 
performance optimization of TEC with space restrictions is a 
very valid assumption. The genetic algorithm is used to 
maximize the cooling capacity and COP of a TEC with space 
restrictions in two different optimization problems. The 
optimization results are validated through finite-element 
simulations using ANSYS®. 


2. DESCRIPTION OF 
COOLER MODEL 


A  THERMOELECTRIC 


The general schematic diagram of a practical single-stage 
thermoelectric cooler is shown in Figure 1 (a). A 
thermoelectric cooler (TEC) consists of many thermoelectric 
(TE) elements. These thermoelectric elements are assembled 
electrically in series. Copper tabs are used to interconnect n- 
type and p-type elements. This array configuration is 
sandwiched between two thermally conducting ceramic 
substrates. Figure 1 (b) is an exploded view diagram of a 
practical TEC system. 

The basic unit of the physical model of a TEC is a 
thermocouple (pair of n-type and p-type semiconductor 
thermoelectric elements). The number of pairs of 
thermoelectric elements may vary from several to hundreds. 
On the one hand, the manufacturing cost of TEC is high, and 
on the other hand, many TE materials are high-priced. Further, 
to predict the performance of a TEC with a heat sink, knowing 
the temperature at important points is quite difficult. Also, the 
thermal resistances in the heat sink, copper conductors, and 
ceramic substrates play a significant role in the total resistance 
to heat flow inthe TEC system. These issues make the 
performance optimization problem challenging to solve. In 
this work, the effects of electrical contact resistance and 
thermal resistance are included. The impact of Joule heat and 
thermal conduction are included as well. 

In this work, to simplify the investigation considering 
thermal resistances, a thermal-resistance model has been 
developed. This model includes thermal resistance of copper 
tabs, ceramic substrates, and cold side heat sink for developing 
a more realistic TEC model. A thermocouple and the 
developed thermal resistance model for this work is shown in 
Figure 2. 
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Figure 2. (a) Thermoelectric couple (b) Thermal resistance 
model 


By applying the electrical analogy of the heat flow to the 
thermal resistance model shown in Figure 2(b), the 
temperatures at the TEC hot surface and the cold surface can 
be expressed as 


[, = Q, (R,,s ot R,, a R.,) + sis (1) 


T= ORR, (2) 
where, 7; and T, are the hot and cold side temperatures (K) of 
n-type and p-type elements. Q, is the heat rejection rate (W) 
from the hot side. Q- is the heat absorption rate at the cold side 
(W), which is referred to as the cooling capacity in common 
usage. To and 7). are the temperatures (K) at the cold surface 
and hot surface of TEC, respectively. Ry; is the thermal 
resistance (°C/W) of the heat sink attached to the hot side of 
TEC, R-, is the thermal resistance (°C/W) of the ceramic 
substrates, and R,, is the thermal resistance (°C/W) of the 
copper tabs. T, is the ambient temperature (K). 

In the current study, some reasonable assumptions are 
considered. 


(a) Heat transfer is assumed to take place along the length 
of TE elements. 

(b) The thermoelectric elements have the same cross- 
section and length. 

(c) Thomson effect is not considered. 

(d) Steady-state condition is prevailing. 


A constant electric current pass through the circuit of 
dissimilar semiconductors. The heat is pumped to one of the 
two sides. It results in making one side cool and another side 
hot. A heat sink attached externally to the hot side ceramic 
substrate dissipates heat to the ambient environment. A 
thermoelectric couple produces cooling or heating effect 
depending on the direction of the electric current. Eq. (3) and 
Eq. (4) shows the heat energy balance at the cold and the hot 
side of the thermoelectric cooler. JT; and T;, correspond to the 
temperature at TE element-copper conductor interface at the 
cold side and hot side, respectively, and used with the same 
reference in each referred equation of this paper. 


kA(T,-T.) 1,,( pL. .r.)| 
=2N| laT oa P +2 
Q. late 1 5 (4 A} (3) 
kA(T,-T.) 1 5( pL. re)] 
=2N| lat. pee py? +2 4 
2, (is is 2 ( A Al (4) 


where, thermoelectric material properties a, p, k are the 
Seebeck coefficient (V/K), electrical resistivity (Qm) and 
thermal conductivity (W/mK), respectively. r, is the electrical 
contact resistance (Qm7). L and A are the length (m) and cross- 
sectional area (m*) of n-type and p-type thermoelectric 
elements, respectively. J is the supplied electric current (A), 
and N is the total number of thermoelectric couples. There are 
three essential terms on the right side of Eq. (3) and Eq. (4). 
The first terms, JaT. and JaT), represent the Peltier heat at the 
cold junction and hot junction, respectively. The second heat 
transfer term kA (T;,—T,)/L is due to thermal conduction. The 
third term % F (pL/A+2r./A) represents the Joule heat 
generation. 

The selection of thermoelectric materials directly affects the 
performance of TEC. The material properties of 
thermoelectric elements are temperature dependent. Bismuth 
telluride (BizTe3) is the popular thermoelectric material used 
in thermoelectric coolers. The material properties of BizTe3 
used in this work are given below, as specified by Fraisse et al. 
[30]. Taveis the average of T, and T),. 
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a=(22224 +930.6T,,, -0.9905T,,, }x10” (5) 
p =(5112+163.4T,,, + 0.6279T,,,)x10-" (6) 
k = (62605 -277.7T,,, +0.4131T.,, )x10“ (7) 


Cooling capacity (Q,) is one of the significant performance 
indexes of TEC, which is used in this study. The Coefficient 
of Performance (COP) is another crucial performance index of 
thermoelectric coolers. Both performance indexes are 
considered in the current study. COP is the ratio of cooling 
capacity to power consumption and defined by the following 
equation. 


Coefficient of Performance,COP = “ 


(8) 


The input electric power (P), as shown in Figure 2(b), can 
be calculated by the following relationship. 


Input Electric Power, P = Q, —Q. (9) 


The cost-competitive and high-performance TEC system 
will pave the way for a promising future of such green devices. 


3. METHODOLOGY 
The various geometrical, material and operational 
parameters affect the cooling performance of the 


thermoelectric cooler. Besides, the restricted maximum area of 
cooling devices, which depends on its application in electronic 
devices, is a_ significant constraint for TEC design. 
Performance optimization is vital to enhance the use of 
thermoelectric coolers in real-world applications. In this study, 
the objective is to maximize the cooling capacity of TEC with 
space restrictions. This paper presents a new approach by 
selecting electric current, length of n-type and p-type TE 
elements and cross-sectional area of TE elements as design 
variables. 


3.1 Optimization of cooling capacity of TEC 


The  single-objective optimization problem for 
maximization of the cooling capacity of TEC is formulated 
mathematically as: 


Maximize Q. 
Subject to 
<J 


(10) 


max 


Tis 
X 


I min 
Lain L 
Anin A < Anna 

Further, the total number of thermoelectric couples (JN) is a 
dependent design variable. Its value depends on the cross- 
sectional area of n-type and p-type thermoelectric elements 


and computed using Eq. (11). 


_ Available area(S )of TEC x packaging density 
2xA 


(11) 


The optimization problem, as mentioned in Eq. (10) has 
been solved using some specific values of parameters. Table 1 
lists the values of the parameters and properties used in this 
work. 


Table 1. Values of parameters and properties 


Description Parameter Value 
Cold surface temperature Tco 293.15 K or 20°C 
Ambient temperature Ta 298.15 K or 25°C 
Heat sink thermal resistance Ris 0.10 °C/W 
Electrical contact resistance re 1x 10° QO m 
Available C.S. area of TEC S 15 mm x 15 mm 
Packaging density PD 80% 
Ceramic thermal conductivity Ker 35.3 W/m°C 
Copper thermal conductivity Kew 386 W/m-°C 


The design variables in the present study are constrained by 
lower and upper bounds. From a practical viewpoint, the range 
for length and cross-sectional area of n-type and p-type TE 


elements is taken as 1.0-2.0 mm and 1.0-2.0 mm’, respectively. 


The range for input electric current is taken as 0.1-3.0 A. The 
dependent design variable N will vary from 45 to 90 as it is 
governed by Eq. (11). The thicknesses of ceramic substrates 
and copper tabs are taken as 0.2 mm and 0.1 mm, respectively. 
The surface area of the ceramic substrate on each side is 
considered identical to the size of TEC. The total surface area 
of the copper tab on each side is considered 90% of the size of 
TEC. Re and Ry are computed as 0.025181 °C/W and 
0.001279345 °C/W, respectively. All these values are taken 
with the help of TEC manufacturing companies’ catalogues. 

Genetic algorithm (GA) is an evolutionary algorithm based 
on natural genetics. The genetic algorithm begins with the 
creation of a population of possible solutions (called 
individuals). Based on the value of the objective function, each 
member of the population is assigned a fitness value. To 
evolve better solutions, new generations are created by 
undergoing selection, recombination, and mutation of 
solutions. The fitness of the new generation is evaluated. This 
cycle is repeated over generations until the stopping criterion 
is met. The objective of GA is to search for an appropriate 
solution for the design problems. This involves maximization 
or minimization of the objective function. 

Genetic algorithm is a population-based optimization 
approach to find optimal or near-optimal solutions. In terms of 
quality and robustness of solutions, GA's capability has been 
widely recognized for providing excellent results on classic 
discrete and continuous optimization problems. The genetic 
algorithm's performance depends on many genetic parameters 
such as population size, crossover, and mutation rate. GA 
parameters play an important role, and a different combination 
of parameters may lead to a significant GA performance 
change. The smaller population size helps faster convergence 
than larger population sizes. The decision on various GA 
parameters and operators are usually selected based on 
recommendations made by GA researchers. 

The real-variable GA employing SBX operator created by 
Deb and Agarwal is used in this study [31]. Table 2 lists the 
values of the GA parameters like population size, crossover, 
mutation & number of generations that are used in the present 
study. The results are reported after multiple runs of GA 
converged to the same best solution. 
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Table 2. Values of GA parameters 


Parameter Value 
Population size 50 
Crossover probability 0.80 
Mutation probability 0.25 
Number of generations 1000 


3.2 Optimization of Coefficient of Performance (COP) of 
TEC 


The objective of the second optimization problem is the 
maximization of the coefficient of performance of TEC. The 
design variables are the same as those selected in the previous 
problem. The fixed values of the parameters and properties are 
identical to the values used in the previous problem and 
described in Table 1. The thicknesses of ceramic substrates 
and copper tabs have the same values of 0.2 mm and 0.1 mm, 
respectively. This new problem is mathematically expressed 
as: 


Maximize COP 
Subject to 
Dinin SL S Lax (12) 
Loin SL S Linax 
A, <A<A 


The goal of this optimization problem is to find the design 
variables within the variable bounds that result in the 
maximum COP of the device. 


3.3 Optimization procedure 


To apply the genetic algorithm to the optimization problems 
described in Eq. (10) and Eq. (12), the fitness evaluation of 
solution vectors is required. However, the procedure for 
evaluating fitness function is slightly tricky for this problem. 
The unknown values of 7; and T, are initially guessed for 
approximately estimate Q. and Q,, The initial guess for T, and 
T. satisfies TEC's prevailing temperature conditions, i.e., 
Tn >Ta and T, <Teo. In principle, these conditions must be 
satisfied. The initial guess will be iteratively modified and 
reach the exact value. Eq. (1) and Eq. (2) are used to calculate 
new values of 7), and 7, that are termed as 7), and T-n. These 
are updated repeatedly to corresponding new values until the 
difference in old values and new values are negligible. Then 
the values of Q. and Q; are accepted. 

A flowchart for GA implementation for these two 
optimization problems are given in Figure 3. 

The brief steps of the fitness evaluation procedure for a 
population individual (solution vector) followed in this work 
are described below. 


(a) 
(b) 


The hot side and cold side temperatures of TE elements 
(T;, and T;) are initially assigned to a guessed value. 
The material properties are estimated using Eq. (5), (6) 
and (7). 

The expected values of Q. and Q, are calculated using 
Eq. (3) and (4). 

Using Eq. (1) and (2) the new values of T;, and T;, are 
calculated. These are termed as Tn and Ten, respectively. 
If the difference of guessed values and new values is 
considerable, then guessed value is updated as 7), = Tin 


(c) 
(d) 
(e) 


and T, = T:n. Go to step (b) and repeat the iteration. 

If the difference of guessed values and new values is 
small, then accept the solution. Take the next individual 
in the GA population to evaluate until all individuals of 
the current generation are evaluated. 


Initialize Population 


Evaluate fitness of 


(f) 


Generation= 
Generation+1 


individuals 


Take first Individual of 
the population 


Guess 7), and T, 


Calculate N from Eq. 
(11) and 
Calculate a,p and k from 
Eq. (5), Eq. 6) and Eq. 
(7) 


Calculate Q, and Q,, 
from Eq. (3) and Eq. (4) 


Calculate Tin and Ten 
from Eq. (1) and Eq. (2) 


Update 
T.=T,» and 
T= Thin 


Is |T4-Thn | < Tolerance? 
and 
Is |7--Ten | < Tolerance? 


Output Values of J, L, and A, 
and N, 0., 01, T, and T;, 


Is the fitness Take next 
ofeach  individu- 


al in the population 
evaluated? 


Individual 
of the 
population 


Apply GA 
Operations: 


Maximum 
Generation 
reached? 


Selection, 
Crossover 
and Mutation 


Report optimal solution and 


fitness vector 


Figure 3. Flowchart for GA implementation 
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4. RESULTS AND DISCUSSIONS 


In the first segment of present work, cooling capacity Q., 
the first performance index of TEC is maximized. The 
algorithm of this study is coded in C language. The GA source 
code is developed by Deb and used in this work [32]. Multiple 
runs of 1000 generations have been repeated, and the best run 
is reported in Table 3 on which algorithm converged several 
times during various runs. 


Table 3. Result of GA based optimization for maximum Q, 


Optimized Optimal Values of Design Variables 
Oc I L A N 
(Dependent) 
8.476807 W 2.993 A 1.0mm 1.607 mm? 56 


At optimal values of design variables, the corresponding 
values of T, and T, are found at 28.59°C and 19.78°C, 
respectively. The hot surface temperature (7j.) of TEC is 
27.84°C. The heat rejection rate (Q;,) at the hot side is 28.401 
W. For the maximized Q,, the value of COP obtained is 0.425. 
It can be observed that L is hitting lower bound while other 
parameters have optimal values without hitting any bound of 
the permitted range. 

To optimize the second performance index of TEC, the 
coefficient of performance (COP) is maximized. The 
boundary conditions and assumptions are similar to those 
considered during the optimization of Q,. This optimization 
problem is solved using the same parameters of GA, as 
mentioned in Table 2. The steps to implement GA in this 
problem are similar to those used in the optimization of 
cooling capacity and shown with the help of a flowchart in 
Figure 3. Several runs of 1000 generations have been 
performed to reach solutions with the highest quality, and the 
best run is reported in Table 4. It is worth mentioning that GA 
converged to the same results in various runs. 


Table 4. Result of GA optimization for maximum COP 


Optimized Optimal Values of Design Variables 
COP I L A N 
(Dependent) 
4.11 0.283 A 2.0mm 1.956mm? 45 


With this maximum COP, the corresponding Q, is obtained 
as 0.745992 W. The corresponding values of 7), and T;, are 
25.11°C and 19.97°C, respectively. The hot surface 
temperature (7),.) of the thermoelectric cooler is 25.09°C. The 
heat rejection rate (Q;,) at the hot side is 0.927 W. The optimal 
values of J and A design variables are unique, while the optimal 
value of L is hitting the upper boundary. It can be seen that 
COP increased significantly, and cooling capacity is just 8.8% 
of Max. Q. obtained, as mentioned in Table 4. It is found that 
a design variable L hits its lower bound for high Q. while for 
high COP, L hits its upper bound. 

From these two results, it is well established that 
maximization of Q, and maximization of COP are obtained at 
a different set of design parameters. Also, maximum Q, does 
not ensure providing optimal COP and vice-versa. This means 
that these objectives are conflicting. The resolution of these 
conflicting design objectives will be Pareto solutions through 
multi-objective optimization if there is no specific objective 
interest. It will be useful to determine a set of solutions that 
will allow the decision-maker to choose among them 


according to the application's requirement. 

The results of this study show that it is possible to improve 
the cooling capacity or COP of the thermoelectric coolers with 
these design variables to be competitive with compressor- 
based cooling devices. The complex impacts of electrical 
contact resistance and thermal resistance deteriorate the TEC 
performance. These factors always need to be included in the 
model for optimization and analysis. 


5. FINITE-ELEMENT SIMULATION 
VALIDATION 


FOR RESULT 


Finite-element simulation is a computational method for 
solving complex engineering problems of the real-world. The 
finite element simulations are performed to validate the 
optimization results of GA. ANSYS® is a useful, common- 
purpose finite-element method tool. It is used to solve a broad 
range of engineering problems numerically. Hence ANSYS® 
is used in the current study. The Thermal-electric module of 
ANSYS® is capable of providing simultaneous solutions of 
thermal and electrical fields. The present work makes use of 
the thermal-electric module for the steady-state analysis of the 
TEC model. A three-dimensional non-linear finite-element 
model is setup. The model in this work is set up with one pair 
of n-type and p-type elements as per the GA result. A new 
approach to incorporate the effect of electric contact resistance 
on the performance of TEC is used in the present study. The 
finite-element simulation includes four additional geometric 
parts termed as ‘Contact’ and used for modelling of the electric 
contact resistance effect. These parts have material properties 
as per the thermo-electric behaviour of electrical contact 
resistance. The contact geometries are created at each end of 
the TE elements. The complete schematic of the TEC model 
for Finite-element simulation to validate GA results is shown 
in Figure 4. 


Contact 


ia 
Current, / 


Figure 4. Schematic of TEC for finite-element simulation to 
validate GA results 


5.1 Finite-element simulation for maximum Q. 


To validate GA predictions for maximum Q,, the length of 
n-type and p-type elements is taken as 1.0 mm, as reported in 
Table 3. The TE elements are of the square cross-section of 
1.27 mm. The distance between n-type and p-type elements is 
0.31 mm. The material properties for the simulation are 
computed at average (Taye) of T, and T, values obtained during 
the GA based optimization of Q,. The finite-element 
simulation input parameters of the modelled TEC are given in 
Table 5. 

To model adiabatic heat transfer from the exposed surfaces 
of TEC, a small convection loss of 0.000001 W/mK was 
applied on all surfaces except the ones on which boundary 
conditions mentioned in Table 5 are specified. The 
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computationally generated mesh, electric voltage, and 
temperature distribution across the finite-element model of the 
thermoelectric cooler are shown in Figure 5. 


Table 5. Finite-element simulation input parameters for 
maximum Q, 


Value per pair 


Description Parameter of TE Elements 
Cooling Capacity Oc 0.1514 W 
Current I 2.993 A 
Temperature (hot side of TEC) Tho 27.84°C 


0 0.001 (m) 
= _—— 
O-OU0S 


(a) 


a ee 
“010613 
0.092864 
0.079598, 
0.066331 
0,053065 
0.039799 
0.026533 
0.013266 
OMin 


0,001 (rm) 
00005 


19.929 Min 


o 


0,001 (rn) 
0.0005 


(c) 


Figure 5. (a) Mesh (b) Voltage distribution (c) Temperature 
distribution in the finite-element model for maximum Qc 


The parameters obtained from finite-element simulation are 
compared with the GA results and reported in Table 6. It is 
observed that the results for a single pair of TE elements from 
GA simulation and those obtained from finite-element 
simulation are in close agreement. The finite-element 
simulation result represents a 3-D solution based on a 
numerical technique, while GA results are based on 1-D 
analytical equations. Hence, the optimization result obtained 
by GA is verified through the solutions of the thermal-electric 
module of ANSYS®. 


Table 6. Comparison of results for maximum Q, 


Parameter GA ANSYS® Remarks 
Ts 20°C 20.25°C es 
On 0.507 W 0.508 W pices fie ae 
P 0.356 W 0.357 W 


5.2 Finite-element simulation for maximum COP 


In this segment, the finite-element simulation for maximum 
COP is performed with ANSYS® software. The steady-state 
TEC model consists of TE elements with 2.0 mm length, as 
reported in Table 4. The TE elements are of the square cross- 
section of 1.4 mm. The distance between n-type and p-type 
elements is 0.38 mm. The temperature-dependent material 
properties are calculated based on the average of 7), and T, 
found during GA based optimization of COP. The input 
parameters of the TEC model for finite-element simulation are 
given in Table 7. 


Table 7. Finite-element simulation input parameters for 
maximum COP 


Description Parameter NaS De pee Ne 
Elements 
Cooling Capacity Oc 0.0165 W 
Current 
T 0.283 A 
Temperature (hot side T 25.09°C 
of TEC) id 


The three-dimensional steady-state TEC model is created, 
and predictions of GA based optimization are tested for 
maximum COP. For this simulation, the mesh, electric voltage, 
and temperature distribution are shown in Figure 6. The finite- 
element simulation results agree well with the GA results. The 
parameters for GA and finite-element simulation results have 
been compared and reported in Table 8. 


o 


0,002 (m) 


0.012227 
0,01048 
0,0087333 
0.0069867 
0.00524 
0,0034933 
0.001 7467 i] 
OMin 


0,002 (m) 
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f=] 


0,002 (rm) 
0,001 


(c) 


Figure 6. (a) Mesh (b) Voltage distribution (c) Temperature 
distribution in the finite-element model for maximum COP 


The ANSYS® result is consistent with GA _ based 
optimization results for the maximization of COP. Hence, the 
optimization result is verified through the solutions of the 
thermal-electric module of ANSYS®. 


Table 8. Comparison of results for maximum COP 


Parameter GA ANSYS® Remarks 
Tco 20°C 19.61°C Value per pair 
Qn 0.021 W 0.021 W of TE elements 
P 0.004 W 0.004 W 
6. CONCLUSIONS 


This paper presents an effective method with a new 
analytical model to improve cooling capacity and coefficient 
of performance of thermoelectric cooler for a specific need. In 
order to analyze more than one factor simultaneously, the 
thermoelectric cooler's current and geometric parameters were 
set to be variables. The described study emphasized to find out 
the optimal values of current, length of n-type and p-type TE 
elements and cross-sectional area of TE elements within size 
restrictions on space. It was found that length, the cross- 
sectional area of thermoelectric elements, and input electric 
current had a great influence on the TEC performance. 
Performance optimizations to maximize cooling capacity and 
to maximize COP were successfully performed by the genetic 
algorithm. The use of this stochastic optimization algorithm 
based on natural genetics theory proved to be the right option. 
The genetic algorithm successfully converged to the same 
optimal results over several runs. The finite-element 
simulations through ANSYS® validated the GA result. 

The work suggests that these design variables should be 
appropriately selected in practical application. Results 
revealed that the relationship between the coefficient of 
performance and cooling capacity is inverse. The maximum 
cooling capacity does not provide optimum COP and vice- 
versa. The smaller length of thermoelectric elements facilitates 
maximum cooling capacity whereas greater length of elements 
obtains maximum coefficient of performance. The best 
performance requires specific values of electric current and 
cross-sectional area of TE elements as per the objective 
requirements. The appropriate optimum results can be 
achieved for any space restriction. This study can guide the 
TEC designers working for some specific cooling targets. The 
use of microprocessor-based control of input power 


parameters to get an optimal cooling with the best possible 
COP under dynamic conditions needs to be explored. 
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NOMENCLATURE 

A cross-sectional area of TE elements, m? 

COP coefficient of performance 

I electric current, A 

k thermal conductivity, W/m K 
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T 


a 


Tave 
Z 
Greek symbols 


or 
p 


length of thermoelectric element, m 
number of thermoelectric couples 

power input, W 

packaging density 

heat absorption rate at the cold side, W 
heat rejection rate from the hot side, W 
thermal resistance of heat sink, °C/W 
thermal resistance of ceramic, °C/W 
thermal resistance of copper, °C/W 
available cross-sectional area of TEC, m7 
temperature at the cold side of elements, K 
temperature at the hot side of elements, K 
temperature at the cold surface of TEC, K 
temperature at the hot surface of TEC, K 
ambient temperature, K 

average of T. and Ty, K 

figure of merit, 1/K 


Seebeck coefficient, V/K 
electrical resistivity, Q m 


